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Black hole solutions in f(R) gravity coupled with non-linear Yang-Mills field 

S. Habib MazharimousavQ and M. Halilsoj|l| 

Department of Physics, Eastern Mediterranean University, 
G. Magusa, north Cyprus, Mersm 10, Turkey. 

It is shown that in the static, spherically symmetric spacetime the problem of metric f[R) gravity 
coupled with non-linear Yang-Mills (YM) field constructed from the Wu-Yang ansatz as source, 
can be solved in all dimensions. By non-linearity it is meant that the YM Lagrangian depends 
arbitrarily on its invariant. A particular form is considered to be in the power-law form with limit 
of the standard YM theory. The formalism admits black hole solutions with single or double horizons 
in which f{R) can be obtained, in general numerically. In 6— dimensional case we obtain an exact 
solution given by f{R) = y/R gravity that couples with the YM field in a consistent manner. 

We dedicate this work to the memory of Yavuz Nutku (1943-2010), the renowned mathematical 
physicist from Turkey. 



I. INTRODUCTION 



With the hopes to explain a number of cosmological problems covering dark energy, accelerated expansion, quantum 
gravity and many related matters, extensions / modifications of general relativity theory gained momentum anew 
during the recent decade. Each extension adds new degrees of freedom and accommodates new parameters apt for 
the sake of better physics. Lovelock gravity JJ, for instance, constitutes one such extension which abides by the 
ghost free combinations of higher order invariants resulting in second order equations alone. Next higher order to 
the Einstein-Hilbert extension in this hierarchy came to be known as the Gauss-Bonnet extension Q which makes 
use of the quadratic invariants. Apart from this hierarchy, arbitrary dependence on the Ricci scalar R which has 
beenpopular in recent times is known as the /(i?) gravity (see e.g. [3| and references therein and for a review paper 
\ see [il). Compared to other theories which employ tensorial invariants this sounds simpler and the fact that the 
V^!} ■ ghosts are eliminated makes it attractive Q. The simplest form, namely f{R) — Ris the well-known Einstein-Hilbert 
' Lagrangian which constitutes the simplest theory of gravity. Given the simplest theory in hand, why to investigate 
ly-j complex versions of it? The idea is to add new degrees of freedom through non-linearities, create curvature sources 
\^ ' that may be counterbalanced by the energy-momentum of some physical sources. Once f{R) — R, passes all the 
. classical experimental tests any function f{R) may be interpreted as a self-similar version of f{R) — R, creating no 
jy-^ ■ serious difficulties at the classical level. Introducing an effective Newtonian constant as a matter of fact plays a crucial 
, role in this matter. At the quantum level, however, problems such as unitarity, renormalizability of the linearized 
theory are of vital importance to be tackled with. 

As an example we refer to the particular form f{R) — R^ {N — rational number) [6]. This admits, among others 
an exact solution which simulates the geometry of a charged object i.e. the Reissner- Nordstrom geometry Q. That is, 
f{R) = R^ behaves geometrically as if we have f{R) = R+ (electrostatic field). Remarkably, the power N plays the 
role of "charge" so that the geometry f{R) = R^ becomes locally isometric to the geometry of Reissner- Nordstrom. 
In a similar manner various combinations of polynomial forms plus ln(l + R), sinR, and other functions of R can 
be considered as potential candidates for f{R). Some of these have already appeared in the literature 7]. Beside the 
case of f{R) = R+ (electrostatic field), and more aptly, cases such as f{R) = R+ (non-minimal scalar field) cases 
also have been investigated [7] . 

Recently, the non-minimal Yang-Mills fields coupled with f{R) gravity has also been studied Q. In this paper we 
show that Yang-Mills (YM) field can be accommodated within the context of metric f{R) gravity as well. To the 
best of our knowledge such a study, especially in higher dimensions which constitutes our main motivation, is absent 
in the literature. While numerical solutions to the problem of black holes in f{R)—YM theory Q started to appear 
in the literature our interest is in finding exact solutions. It should be added also that the class of black holes in 
f{R) gravity can be distinct from the well-known classes such as Myers-Perry (in 0). Some classes of black holes 
in this paper also obey this rule since they are not asymptotically flat in the usual sense. We show that in most of 
our solutions asymptotically (i.e. r — > oo) an effective cosmological constant can be identified which depends on the 
dimension of spacetime, YM charge Q and the integration parameters. Herein, we do not propose an f{R) Lagrangian 
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a priori, instead we determine the f{R) function in accordance with the YM sources. By using the Wu-Yang ansatz 
for the YM field [HI it is shown that a general class of solutions can be obtained in the f{R) gravity where 
the geometric source matches with the energy-momentum of the YM field. Let us note that the Wu-Yang ansatz 
works miraculously in all higher dimensions which renders possible to solve f{R) gravity not only in d = 4, but in all 
d > 4 as well. Further, the zero trace condition for the energy-momentum is imposed to obtain conformal invariant 
solutions which constitutes a particular class. Although our starting point is the non-linear YM field in which the 
Lagrangian is an arbitrary function of the YM invariant our main concern is the special limit, namely, the linear (or 
standard) YM theory. Power-law type non-linearity is a well-known class which we consider as an example and as 
we have shown elsewhere the choice of the power plays a crucial role in the satisfaction of the energy conditions, 
i.e. Weak, Strong or Dominant. Essentially, conformal invariant property is one of the reasons that we consider 
Power- YM (PYM) field in higher dimensions. The implementation of zero trace condition for the energy-momentum 
tensor becomes relatively simpler in this PYM class. As in the example of Born-Infeld electrodynamics case which 
plays crucial role in resolution of point like singularities, in analogy, similar expectations can be associated with the 
non-linear version of the YM theory. Such non-linearities resemble the self-interacting scalar fields which serve to 
define different vacua in quantum field theory. More to that, in general relativity the non-linear terms effect black 
hole formation significantly, it is therefore tempting to take such combinations seriously. It is our belief that with the 
non-linear YM field we can establish effective cosmological parameters to contribute, in accordance with the energy 
conditions cited above, to the distinction between the phantom and quintessence data of our universe. This is a 
separate problem of utmost importance that should be considered separately. We show that, the choice f{R) — VR 
in 6— dimensions yields an exact solution for f{R) gravity coupled with YM fields which is non-asymptotically flat / 
non-de Sitter in the sense that it contains deficit angles at r — )■ oo. Other classes of solutions that are asymptotically 
de Sitter, unfortunately can't be expressed in a closed form as f{R). 

Organization of the paper is as follows. In Section II we introduce our theory of non-linear YM field coupled to f{R) 
gravity and give exact solutions. Section III specifies the non-linearity of YM field to PYM case in all dimensions. The 
First Law of thermodynamics in our formalism is discussed briefly in Sec. IV. We complete the paper with Conclusion 
which appears in Sec. V. 

II. NON-LINEAR YM FIELD IN f{R) GRAVITY 
We start with an action given by 



S= d' 



(1) 



in which / (i?) is a real function of Ricci scalar i?, L (F) is the non-linear YM Lagrangian with F = jtr 
The particular choice L (F) — —j^F will reduce to the case of standard YM theory. Here 

F'-"^ ^^Fj^l^dx^" Adx" (2) 

is the YM fleld 2— form with the internal index (a) running over the degrees of freedom of the YM non-abelian gauge 
fleld. Our unit convention is chosen such that c = G = 1 so that k — 8tt. 

Variation of the action with respect to the metric gives the field equations as 

fnR; + (ofji - i/) 5"^ - V^V.h - (3) 
where fj^ ^ ^ and Ofn - V^V^/^ = i 9^ (^3^9^) Further, V'WJr = 5""^ (/fl).,.. = 



ifR)u 



The trace of the field equation implies 



fRR+{d-l)afR-^f = KT (4) 



2" 



in which T = T!'. 



The energy momentum tensor is chosen to be 



T; = L [F) S; - tr (^f(-Jf(^>'^) Lf (F) (5) 



in which Lf (F) - 



dF 



The YM ansatz, following the higher dimensional extension of Wu-Yang ansatz, is given by 



A(") = ^C^f^Q.) x'dx^, g = YM magnetic charge, = ^a;^, 



2 < j + 1 < « < - 1, and 1 < a < (d - 2) (d - 1) /2, 
xi — rcos6d-3sm6d-4:...sin0i, X2 = r sin^^-s sin0(i-4--- sin^i, 
x^ — r cos 6*^-4 sin ^d-s--- sin 01, a;4 = rsin6'd-4sin0£j_5...sin0i, 

Xd-2 = r cos 6*1, 

in which C'jflj^jc) ^® non-zero structure constants [loj . The spherically symmetric metric is written as 



ds^ = - A [r) dr + — — + r^dni^2, 



A{r) 



where A (r) is the only unknown function of r and 



d-2i-l 
2 , V- TT „:„2 /). ,q2 



dni_2 = rf^t + E n Sin' 0, doi 

i=2j=l 



with 



The YM equations take the form 



< Od-2 < 27r, < e*,, < TT, 1 < i < d - 3. 



*F(°)Li. (F) + -C\^}^^^^Lf (F) AC") a* F('=) = 0, 



where * means duality. For our future use we add also that 



4 \ f^i^ 



4^4 



and 



while 



(d-3)g2 



From (5) the non-zero energy momentum tensor components are 



n = l = t;, 



The trace of Eq. (5) implies 



which yields from Eq. (3) 



{d~i)Q\ 
Iff, - L -. Lp- 



T = d.L- AFLf 



!=-, [IrR + (rf - 1) Dfn - K {d.L - AFLf) 
d 



To write the exact form of the field equations we need the general form of Ricci scalar and Ricci tensor which are 
given by 

„ ^ r'A" + 2(d-2)rA' + (d-2)(d-3){A-l} ^ ^^^^ 
in which a prime denotes derivative with respect to r. Overall, the field equations read now 



2 r y V 2 

fn ( + + (af^_K^_ v-v./h = (21) 



2 r J \ 2' 

rA' + {d -3){A-l) \ ( _1 



Herein 



□/k = -^a, (7^5'^) /fl A'/;^ + Af^ + ^^^/ii, (23) 

V*V = = 5" {fR,t,t - TTtfR,m) = ^^'/fl, (24) 

V'V./ij = = 5'" - r™/K,„) = Af^ + 1^7^, (25) 



and 



V^^V.Jfl = = 5'*'* [fRAfi. - ^'SlejR.m) = ^fk- (26) 

The tt and rr components of the field equations imply 

VV./fl = V*Vt/H (27) 

or equivalently 

f'A = 0. (28) 

This leads to the solution 

fR = ^ + Vr (29) 

where ^ and rj are two integration constants and to avoid any non-physical case we assume that 7],^, > 0. The other 
field equations become 

IrA" + (d-2)A'\ 1 (rf-2)^ 1 



In a similar manner the 6i9i and components yield 

- r'^A" - 

~2r^ ) ' \r 2)'' '"{(1-2) 



, , 2(d-3)(A-l)-r2A"-(d-4)rA'\ [A A'\ 4 
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Equation (16) can equivalently be expressed by 

2 ' 



JrR +{d-l)[ A't] + ^ — "^Arj ] - K {d.L - AFLp) 



(33) 



Here we comment that in the Hmit of hnear Einstein- YM (EYM) theory one may set L = —j^F, Lp = ^j^- ^ = 0, 
and ^ = 1 to get //j = 1 or equivalently f = R and consequently 

2 (d - 3) (A - 1) - r^A" -id -A) rA' = -Ifc^I^. (34) 

This admits a solution as 

which was reported before Here we note that m is an integration constant related to mass of the black hole. 

III. PYM FIELD COUPLED TO / (R) GRAVITY 

A. General integral for the PYM field in f{R) gravity 

Our first approach to the solution of the field equations, concerns the PYM theory which is a particular non-linearity 
given by the Lagrangian L = —j^F'^, in which s is a real parameter [l^. The EYM limit is obtained by setting s = 1. 
The metric function, then, reads as (^ = 0) 



d-3 I A„2 m {d-l)(d-2)'^2^ {d~3)''i'' Q^lnr „ _ d-1 

A{r)^{ "-'^ 2^^d ' ' ^ (36) 

±:1 ^ \r-'^ ™ 4''-'s(d-2)°-^(d-3)°Q"° / 

d-2 ' r<i-2 (4s+l)(d-4s-l)?)r«'s-i ' * 7^ 4 

in which A and m arise naturally as integration constants. Obviously A is identified as the cosmological constant 
while m is related to the mass. The fact that our metric is asymptotically de Sitter seems to be manifest only with 
deficit angles at r — oo. We add that in order to have an exact solution we had to set ^ = 0, which means that in 
this case we were unable to obtain the f{R) = R gravity from the general solution. 
Using the metric solution we also find / (R (r)) and R (r) as 

d-l ei-5 d-1 

n„d-3 (d-3) — (rf-2)— Q— _ d-l 

f{R{r)) = { ' ' 4 (37) 



9„d-3 (d-3Y(d-2r-\is-d+2) « ^ d-l 
JJ^l 733", S ^ ^— 



and 



d-l d-5 d-l 



^-Ad(d-i)- '"-^^('^-^]T^'~^^^^T , s = ^ 

^ - Ad (d - 1) - '^-'^+1-^1)'"^°'^ ;Cr, . ^ ^ 
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We recall from equation (29) that 



or equivalently 

djjdr 
dR/dr 



rjr. (40) 



As it is seen from the expressions of R (r) and / (r) it is not possible to eliminate r to have the exact form of / [R) , 
instead we have a parametric form for / (i?) . 
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Among all possible cases, we are interested in the condition 4s — d + 2 = 0. Since this particular choice brings 
significant simplifications in (37) and (38). Table 1 shows for which values of s and d this is satisfied. 



d = 


5 


6 


7 


8 


9 


10 


d 


s = 


3 
4 


1 


5 
4 


3 
2 


7 
4 


2 


d-2 
4 



It is not difficult to observe that with these specific choices, plus A = 0, we obtain 
in which the constant fi^ is defined by 



/io = 2r]\/d- 3. 
Accordingly the metric function A (r) takes the form 

d — 3 m 



A{r) 



d-2 



d-2 d—2 j_9 

{d-2)— {d-3)—Q^ 
{d - 1) 



with the scalar curvature 



R{r) 



d-3 



Since the constant term 7^ 1, in (43) our solution for r — !• cx) is given by 



ds^ + ^,.2 ^ (^-)2 ^^2__^ 



(Table 1) 



(41) 



(42) 



(43) 



(44) 



(45) 



where t — ^/at, r — ar, (a — ^35). This may be interpreted as a deficit angle at r — >■ 00. It can also be seen easily 
from Table 1 that for the linear YM theory (s = 1) in c? = 6, with ?y = f{R) — VR yields an exact solution. 



B. Thermodynamics of the black hole solution 

The black hole solution given by (36) admits horizon(s) provided 

A(r^) = 0, 

which implies 

d-3 d-2 . d (d-l)(d-2)V-QV-(d-3) V- 1 o - 

d^^h + ^^^h d-5 inr/i, s — — 

d-2 , . d _ i^-''s(d~2y~\d-3yQ^'' I 
d-2'h -r^^'h (4s+l)(<i-4s-l),,r*=-'*+i' 



d-1 
4 



The standard definition of Hawking temperature 



yields 



ad(d-3)^/"((rgA+l)d-3)-V2(d-2)^(d-l)r^-''QV(^)4r;, ^ ^ d_l 

Th ^ { (d_3)i^rfc,,d ' 

[,,(4s+l)(d-2)(r''Ad+d-3)-4(ii^)°(d-3)''Q^'^r^-^°] 

47rr;(d-2)(4s+l)rh ' 7^ 4 

It is known that the area formula S = in /(i?) gravity becomes [l3 - [l6| 

a _ f I 

4G •'^''■='■'1 



(46) 



(47) 



(48) 



(49) 



(50) 
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in which 



(51) 



and 



Ah = 



d-l 



(52) 



where is the radius of the event horizon or cosmological horizon of the black hole. Using S with the definition of 
the heat capacity in constant charge we get 



Co = T, 



H 



( dS 



(53) 



2s„-4s+l 



^^7? jd - if r"-^ r,id~2){s + 1) {Ar^d + d ~ 3) - As jd ~ 3)' Q^'r 

AT (^) rj{d-2){s + i) {Ar^d - d + 3) + 16s^ [d - 3)' Q2.^-4s+i ' 

A through analysis of the zeros / infinities of this function reveal about local thermodynamic stability / phase 
transitions, which will be ignored here. 



C. A general approach with s = 1 

1. d> 6 



In this section, for the linear YM theory (s = 1), we let ^ to get non-zero value and attempt to find the general 
solution. As one may notice the case of o? = 5 is distinct so that we shall find a separate solution for it but for d> 6 
the most general solution reads 



A{r) 



1 



(d - 3) ( 



{d - 5) 



{-lfr^-f'^-\d~l)m\xi 



+ Pd-7 (7) , 



in which 7 = | with the following abbreviations 



Ai 
A2 
A3 
A4 



d-l 



A{d~3) Q2 
3(d-5)^V ' d-2^r 



2 r 



rf-4 



In 



{d~2) {d~b) e 

4(d-3)g2 d-l m 
-r - 



d-2C 

d-l m 2 (d - 3) Q2 

' d-3(\d-5 ^ (d-b) ~^ 



rf-5 ^4 d-Afr'i-e' 



4(d-3)Q2 

— ; — ; In 

d-5 



d-l I 



d — 5 r 



d-7 ■ 



(7) - (-l)''"'(d-l)m7 



rf— 6 tcrm(s) 



d-& 



(54) 



(55) 



d>7. 



Using these results, we plot Fig. 1 which displays (for d — 6), A{r) and f{R) for different values of rj. 
It can easily be seen that in the limit 77 — > and ^ = 1, the metric function reduces to 



A{r) = l 



{d-3)Q^ 
{d - 5) r2 ' 



(56) 



which is nothing but the well-known black hole solution in f{R) — R, EYM theory [ll|. On the other hand for 
ry 7^ 7^ ^, it is observed that asymptotical fiatness does not hold. 
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To complete our solution we find the asymptotic behavior of the metric function A{r). As one observes from (54), 

at r — oo A{r) becomes 



A(r)~l + ^r2, 



(57) 



in which 



A 



eff 



(d-2) (d~5) 



4(d-3) 
d~5 



QV + (-l)S''-'(c^-l)e 



m 



7^ In . 



(58) 



As one may see in Fig. 1, we add here that j] plays a crucial role in making the metric function asymptotically de 
Sitter, anti - de Sitter and asymptotically flat {\r]\ < 1, I77I > 1 and |?7| = 1 respectively). 

a. Thermodynamics of the BH solution in 6— dimensions In this part we would like to study the thermodynamics 
of the solution (54) and compare the result with the case of linear gravity f{R) = R. As one can see from the form of 
the solution (54) , wc; arc; not able to study analytically in any arbitrary dimensions d, and therefore we only consider 
d = 6. The metric solution in d = 6 dimensions is given by 



Mr) 



^ _ m _ 3Q2 



1 



40^ 



+ 



C + w 


5 m 


r 


"2^ 


5 m^ 




+ 2^y 





5 TO 1 r 
~ 2| 

e 



+ 



(59) 



r^ln 





- 5r^ 


r 





and therefore the Hawking temperature is given by 
3 3 Q2 ^^ / ^ 3 1 



+ 



and the specific heat capacity reads 



yiQirrl 16 tt/^ V 



f U'SQ^ 120r,,lnr;, + 139r;, \ 
GAwrh 19277 / 



Cq = 



(3r2-9Q2) 



8 7rV|(3r^-17QM + 7Q4) 

9 {rl-3Q^f 



(60) 



(61) 



First we comment herein that, to get the above result we considered ^ = 1. Second wc add that, by 77 = wc; get the 
case of EYM black hole in i?— gravity. In the case of pure i?— gravity we put 77 = and Q = which leads to 



, and Cq = - 



(62) 



which are the Hawking temperature and Heat capacity of the 6— dimensional Schwarzschild black hole. Divergence in 
the Heat capacity for particular YM charge and therefore a thermodynamic instability is evident from this expression. 



2. d = 5 



As we stated before, dimension d = 5 behaves different from the other dimensions. The metric function is given in 
this case by 



in which 



Ai = 
A2 = 
A3 = 
A4 = 



9fr 
1 

3^ 
2r 
74 



TO 

A{r) = 1- ^ - In'' + ??Ai + rj'^ A2 + ry^Ag + r]^Ai 



[24Q2 In r + 12to + 2Q^ - 6^r 



2er2ln-^ + '^'" io^2,_ ,^2 



r 

[2to^ - 3Q2 + 4g2 In 



- 12g2lnr-3Q^-6TOC 



2r2 



4Q^ In r In + 4Q2dilog 



^ + ^^V(2TO^-Q^)ln^^ + ^'' 



(63) 



(64) 
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Herein m is an integration constant and 



dilog (a;) 



Int 
1-f 



-dt 



(65) 



is the dilogarithm function. Here also the EYM hmit with — > and ^ = 1 is obvious. 

Similar to the higher than 6— dimensional case we give here also the asymptotic behavior of the metric solution (65) 
as r — > oo, 



(66) 



in which 



A 



3mf In r;-Q^ In ( ^ ) In (C?/) 



ln77 + 27r2Q2 



(67) 



is the effective cosmological constant. 



D. Black holes with a conformally invariant YM source 

One of the interesting choice for s in Einstein-Power-Maxwell theory - which has been considered first by Hassaine 
and Martinez [l^l - is given by s = | (for all d > 4) which is conformally invariant. This choice yields a zero trace 
for the energy momentum tensor in any dimensions, i.e., T = Tj^ = 0. In EPYM case also s — j leads to a traceless 
energy momentum tensor and a metric solution for the field equations with arbitrary values of f and t] is given by 



d-2 



In 



e(l+7r) 



27 
d-2 



r+g(7) (68) 



in which 



d-2 



q{j)^i-lfid-l)mY-'J2 



i-irr 



k=l 



7 



k-ll 



(69) 



and 7 = |. Fig. 2 (d = 5, s = |) and Fig. 3 (d = 6, s = |) depict A{r), R{r) and f{R) which relate the conformally 
) cases for different rj values. For r — >■ oo, it can be seen easily from Eq. (68) that we have an effective 



invariant (s = 
cosmological constant term, given by 

Ae// = 3 



+{-lfid-l)mY-^ 



d-2 



ln|77| 



(70) 



We note that as a limit, once rj ^ (or equivalently 7 -> 0) and ^ — > 1 the solution reduces to 



A{r) 



m 



-4 {d-2) 



d-3 



„d-2 



(71) 



which is the metric function in Einstein-PYM theory in /(i?) = R gravity. Determination of horizons and thermody- 
namical properties in this limit is much more feasible in comparison with the intricate expression (68). To complete 
this section we give the Hawking temperature and specific heat capacity for d — 5 which read 



^ 2rl-i/2AQ^ ^ 2^240^- 17rl + 10^/2AQ^ + Urjlnr,, ^ ^ 



(72) 



and 



Cr 



24Q10 _ ) Sn'^r'l (4 {/2AQ^rl - - 2^/&Q' 



2y2AQ 



10 



(2^/2AQ 



10 _ ^3 
' h 



(73) 
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1. Constant d— dimensional curvature R = Rq 

G. Cognola, et al in Ref. [l^ have considered the constant four-dimensional curvature R — Rq in pure / (R) 
gravity, which imphes a de-Sitter universe. Here we wish to fohow the same procedure in higher dimensions in / (i?) 
gravity coupled with the non-minimall PYM field. As stated before, in order to have a traceless energy momentum in 
d— dimensions we need to consider the case of conformally invariant YM source which is given by s = I in the PYM 
source. In 4— dimensions s = 1 is satisfied automatically for the zero trace condition. 

We start with the trace of the Eq. (4) which leads to 



nRo) = i^f{Ro) (74) 



and therefore the field equations (3) become 



G'; + A,ffS;^^f; (75) 
with the effective cosmological constant and energy momentum tensor as 

(d - 2) Rq 2Ro 

^ 2d ' ^^ = jmd^^- 

Now, we follow [H, [iB] to give the form of the entropy akin to the possible BH solution. As we indicated in Eq. (50) 
the entropy of the modified gravity with constant curvature is given by 

S-^h. (77) 



which after considering (78) it becomes 



5 = ^/(i?o). (78) 



Since our main concern in this paper is not the particular class oi R = Rq =constant curvature space time we shall 
not extend our discussion here any further. 



IV. FIRST LAW OF THERMODYNAMICS 



In this section we follow Ref. [l7j to find a higher dimensional form of the Misner-Sharp energy [18] inside the 
horizon of the static spherically symmetric black hole in f{R) gravity. The corresponding metric is given by (7) and 
the horizon is found from A (rh) = 0. The field equations (3) may be written as 



(79) 



where is the Einstein tensor and is a stress-energy tensor for the effective curvature which reads 



1 



1 



(80) 



At the horizon tt and rr parts of (79) imply 



d-2 {d-2) id-3) 
^^^^ 



fB.^^[ + ^ [(/ - RfB) - A'f 



1 1) ■ 



(81) 



which upon multiplying by an infinitesimal displacement dr^ on both sides can be reexpressed in the form 



An 



2TiAh 



{d-2) (d-3) 



Ir + (/ - RIr) 



Ahdrh = AhT^drh- 



(82) 
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We add here that all functions are calculated at the horizon, for instance A' — i^^lili 
suggests that 



dr 



dE^ — 
2k 



{d-2) (d-3) 



Ib. + (/ - RfR) 



Ahdrh 



The latter equation 



(83) 



in which E is the Misner-Sharp energy in our case. Therefore (82) becomes 

TdS -dE = PdV 

A' 



(84) 



where we set Hawking temperature T — j^, entropy of the black hole S = Jr, radial pressure of matter fields at 
the horizon P = T^ = T^ and finally the change of volume of the black hole at the horizon is given by dV — Ahdrh- 
The exact form of the Misner-Sharp energy stored inside the horizon may be found as 



E=^ 
2k 



(d-2)(d-3) 



fR + (/ - RIr) 



Ahdrh 



(85) 



in which the integration constant is set to zero (to read more see Ref.s daiil). 

As an example we study the case of PYM field in f{R) gravity in Sec. III-A. Also to have an exact form for f{R) 
we employ the metric (43) which corresponds to f{R) ~ a^qV^- Eq. (82) yields, 



f2TTAh \ 




[ K '^'^V 






2k 



(d-2)(d-3)^^ ??(rf-3) 



rh 



rh 



Ahdrh = Ah 



-1 
47r 



{d-2) (d- 3)( 



drh 



(86) 



in which R 



d~3 



has been used. Now, this equation leads to 



A' 



(d-3) 



1 



rh 



vid-1) 



{d-2) {d-3)Q^ 



By taking derivative of (43) and substituting for m in terms of rh the foregoing equation easily follows. 
Finally, one can see that the total energy is expressed by 



E 



r]{d-3) {d-1) 
2k {d - 2) 



Ah- 



(87) 



(88) 



V. CONCLUSION 



An arbitrary dependence on the Ricci scalar in the form of f{R) as Lagrangian yields naturally an arbitrary 
geometrical curvature. The challenge is to find a suitable energy-momentum that will match this curvature by solving 
the highly non- linear set of equations. For a number of reasons it has been suggested that f{R) gravity may solve 
the long-standing problems such as, accelerated expansion and dark energy problems of cosmology. Richer theoretical 
structure naturally provides more parameters to fit recent observational data. We have shown that in analogy with 
the electromagnetic (both linear and non-linear) field, the Yang-Mills field also can be employed and solved within the 
context of f{R) gravity. So far, f{R) as a modified theory of gravity has been considered mainly in d = 4, whereas 
we have been able in the presence of YM fields to extend it to d > 4. In addition to the parameters of the theory 
the dimension of space time also contribute asymptotically to the effective cosmological constant created in f{R) 
gravity. Admittedly, out of the general numerical solution technically it is not possible to invert scalar curvature R{r) 
as r = r{R) and obtain f{R) in a closed form. This happens only in very special cases. In particular dimensions 
and non-linearities we obtained black holes with single / multi horizons. From the obtained solutions for PYM field 
coupled f{R) gravity we can discriminate three broad classes as follows: 

i) the asymptotically flat class in which 77 = 0, ^ = 1. This class was already known [ll| . 

n) the asymptotically de-Sitter / anti-de Sitter class corresponding to 77 7^ 0, f = 1 (s = 1). 

iii) the non-asymptotically flat / non-asymptotically de Sitter class for 77 ^ 7^ ^, s = |. 

Our solutions admit black hole solutions with single / multi horizons. In the proper limits we recover all the 
well-known metrics to date. The case (ii) at large distance limit exhibits deflcit angle as shown in Eq. (45). 

Conformally invariant class with zero trace of the energy-momentum tensor, is obtained with the PYM Lagrangian 
L {F) = —j^{F)^^^ in d = 5. In general, the power of F becomes meaningful within the context of energy conditions 
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and causality. By introducing effective pressure Peff and energy density p^f^ througli Pg// = ^Peff ^'^'^ using tlie 
PYM fields in energy conditions lu factor (i.e. whether cj < — 1, or cj > — 1) can be determined as a cosmological factor 
[3]. This will be our next project in this line of study. It may happen that, certain set of powers eliminate non- physical 
fields such as phantoms and alikes. As far as exact solutions are concerned a remarkable solution is obtained in the 
case of standard YM Lagrangian L (F) — —-^F with c? = 6 in f{R) = VR gravity which automatically restricts the 
curvature to i? > 0. 
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Captions: 

Fig. 1: The plot of 6— dimensional f{R) (Fig. la) and A{r) (Fig. lb). We choose A = 0, and four different values 
of 77 {tIatVetVc ^ii'i Vd) depicted as plots A, B, C and D. From Fig. lb it can be seen that in A, B and C we 
have single, while in D double horizons. 

Fig. 2: The 5— dimensional plots of A{r), f{r) and R{r) from Eq. (68), for a variety of parameters given in Fig.s 
(2a-2d). Since s = | in this particular case, the source is the PYM field with Lagrangian L - These are all 

black hole solutions with inner and outer horizons. A general analytic expression for /(i?) seems out of our reach. 

Fig. 3: The plot of the metric function A{r) corresponding to the conformally invariant case from Eq. (68) in 
d = 6 and s = |, for a set of rj parameters. Black hole formations with single / double horizons are explicitly seen. 
Specifically, Fig. 3a for < < 0.80 and Fig. 3b for 0.9 < 77 < 1.0. 

Table 1: The table for d versus s that satisfies the condition 4s — c? + 2 = 0. The reason for making this choice is 
technical for it simplifies the expressions in (37) and (38) to great extend. 
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